Dynamical quantum phase transitions in the dissipative Lipkin-Meshkov-Glick model 

with proposed realization in optical cavity QED 



S. Morrison 1,2 and A. S. Parkins 3 

institute for Theoretical Physics, University of Innsbruck, Innsbruck A-6020, Austria 
2 Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, A-6020 Innsbruck, Austria 
Department of Physics, University of Auckland, Private Bag 92019, Auckland, New Zealand 

(Dated: February 4, 2008) 

We present an optical cavity QED configuration that is described by a dissipative version of the 
Lipkin-Meshkov-Glick model of an infinitely coordinated spin system. This open quantum system 
exhibits both first- and second-order non-equilibrium quantum phase transitions as a single, effective 
field parameter is varied. Light emitted from the cavity offers measurable signatures of the critical 
behavior, including that of the spin-spin entanglement. 

PACS numbers: 42.50.Fx, 42.50.Pq, 03.65.Ud, 73.43.Nq 



Remarkable advances with trapped, ultra-cold atomic 
gases have opened up exciting new avenues of research 
into strongly interacting many-body quantum systems 
PP. Exquisite control of both motional and electronic de- 
grees of freedom of cold atoms can enable one to "tailor" 
atom-atom interactions and thereby implement a vari- 
ety of systems that exhibit, in particular, quantum crit- 
ical phenomena, i.e., transitions between distinct quan- 
tum phases, driven by quantum fluctuations, in response 
to variations of an effective field or interaction strength 
around some critical value. 

Recently, important insights into such transitions have 
been obtained from theoretical studies of the quantum 
entanglement properties of critical spin systems (see, e.g., 
HHHSISHISlGSlin])- Bipartite entanglement mea- 
sures characterizing entanglement between a pair of spins 
(e.g., the concurrence) or between two blocks of spins 
(e.g., the entanglement entropy) can display marked crit- 
ical behavior and scaling at quantum critical points. In 
this context, a simple but very useful example is the 
Lipkin-Meshkov-Glick (LMG) model [TT], which is de- 
scribed by the Hamiltonian 

Hlmg = -2hJ, - (2X/N)(Jl + 7 J 2 ), (1) 

where {J x , Jy, Jz} are collective angular momentum op- 
erators for N spin-1/2 particles, h and A are effective 
magnetic field and spin-spin interaction strengths, re- 
spectively, and 7 G [— 1,1] is an anisotropy parameter. 
This system, in which each spin interacts identically with 
every other spin, exhibits critical behavior at zero tem- 
perature; in particular, either first- or second-order equi- 
librium quantum phase transitions may occur, depending 
on the choice of A and 7, as the ratio h/Xis varied across 
a critical value [BJ. Notably, the second-order transition 
involves a change from a unique ground state (normal 
phase) to a pair of macroscopically displaced degener- 
ate ground states (broken phase). Entanglement in the 
system displays the above-mentioned critical behavior, 
reaching, in particular, a pronounced maximum at the 
critical point [51[BJ|7]- 



Given these interesting and topical features of the 
LMG model, it follows that the physical realization of a 
system described by such a model would provide a valu- 
able test bed for studies of quantum critical phenomena 
and entanglement. Here we propose an open-system (i.e., 
dissipative) version of the LMG model based on the col- 
lective interaction of an ensemble of atoms with laser 
fields and field modes of a high-finesse optical resonator. 
In the spirit of a recent proposal for realizing the Dicke 
model [12], our scheme employs Raman transitions be- 
tween a pair of atomic ground states and the relevant 
energy scales (e.g., h, A) are set by light shifts of the 
atomic levels and Raman transition rates and detunings. 

The open nature of this system, a consequence of the 
external driving fields and cavity mode losses, introduces 
a number of important differences from, and, arguably, 
advantages over, the closed, Hamiltonian LMG system: 
(i) thermal equilibrium phase transitions are replaced 
by dynamical, non-equilibrium phase transitions, (ii) the 
cavity output fields offer quantitative measures of proper- 
ties of the collective-spin system, including entanglement, 
in the critical regime, and (iii) it is possible to observe 
both first- and second-order quantum phase transitions 
as a single effective field parameter, h, is varied. 

We consider N atoms coupled via electric dipole tran- 
sitions to three laser fields and to a pair of independent 
(e.g., orthogonally-polarized) optical cavity modes. The 
atomic level and excitation scheme is shown in Fig. [I] to- 
gether with a possible ring-cavity setup. At the location 
of the atoms, the cavity and laser fields are copropagat- 
ing traveling waves, with sufficiently broad beam waists 
so as to ensure homogeneous atom- field couplings. These 
fields combine to drive Raman transitions between two 
stable electronic ground states of the atoms, |0) and |1) 
(energies too — and uj±, respectively, with h — 1) via 
the excited atomic states \r) and \s) (energies ui r and 
lu s ). The laser fields have optical frequencies u> r o, w s o, 
and u) r i, and couple to the atomic transitions with Rabi 
frequencies fl r o, f^o, and f2 r i. Cavity field a, at fre- 
quency uj a , couples to the transitions |0) «-> |r) and 
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|1) <-> |s) with strengths g r o and <? s i, respectively, while 
cavity field b, at frequency uj b , couples to the transitions 
|0) <-> |s) and |1) «-> |r) with strengths <? s o and g> r i, re- 
spectively. As drawn in Fig. [TJ the level scheme would 
apply e.g., to 6 Li, with the ground magnetic substates 
|F = 1/2, m = ±1/2) as |0) and |1) and a magnetic field 
perpendicular to the cavity axis to provide a splitting 
lui between these states. Modes a and b would be or- 
thogonal, linearly-polarized cavity modes, with mode a 
polarized along the direction of the magnetic field. 



with J x — ( J + + J_ ) /2 and 




f^rO ,^. s 



FIG. 1: (a) Atomic level and excitation scheme, (b) Poten- 
tial ring-cavity setup. The laser fields (dashed lines) are at 
frequencies that are not supported by the resonator, but can 
be injected through one of the resonator mirrors so as to be 
co-propagating with the cavity fields through the ensemble. 

The atom-light detunings A r = uj r — lo t .q and A s = 
lu s —lo s q are taken to be much larger than any dipole cou- 
pling strengths, atomic linewidths, or cavity loss rates. 
This enables us to adiabatically eliminate the states |r) 
and \s) from the dynamics and neglect the effects of 
atomic spontaneous emission. Additionally, as depicted 
in Fig. [TJ we assume that only three distinct Raman 
transitions are of significance (i.e., resonant or roughly 
resonant); i.e., we retain only those Raman processes 
that cause a change in the electronic state of the atoms 
(|0) — > |1) or |1) — > |0)) and also involve transfer of a 
photon from a laser field to a cavity mode or vice-versa. 
All other possible Raman processes are assumed to be far 
off-resonant and therefore negligible. Finally, taking the 
wave numbers of the laser and cavity fields to be essen- 
tially equal, and introducing the collective spin operators 

J. = iEjLidiiXiil-lOiXOjl), J+ = Ef=i|i,->vU 

and J_ = ( J+y , we can derive a master equation for the 
cavity modes and ground-state atoms in the form 



Pe = -i[H g ,pg\ + K a D[a]p g + K b D[b]p g , 



(2) 



where D[A}p = 2ApA* -A* Ap-pA^A, 

n a ,b are the cavity 
field decay rates, and (omitting constant energy terms) 

H s = u Jz + S a a*a + 5 b tfb+25~J z a}a + 28 b ~ J z tfb 
+ ^J x (a + a^) + ^L(J_b + J + tf), (3) 
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where = (lu s q — Lo r \)/2 ~ u>x, and we have assumed 
the two Raman transitions involving mode a to occur at 
the same rate A Q . 

We now assume (k 2 + Sf ) 1 ^ 2 » A a , A b , lu . In this limit, 
the cavity modes are only weakly or virtually excited and 
may also be adiabatically eliminated to yield the follow- 
ing master equation for the reduced density operator, p, 
of the collective atomic system alone: 



r 



N 



jD[2J x 



^D[J + ]p, (5) 



with h = -w /2, A = 2X 2 a S a /( K 2 a + 8%), and T t = 
\fni/(K 2 + 8 2 ) (i — a,b). Note that in deriving ^ we 
have also assumed that K b ^S> 8 b ~ 0. If we then take 
<5 a n a and T a <C T b , then the role played by each 
cavity mode in relation to the atomic system is quite dis- 
tinct. Specifically, mode a mediates the collective spin- 
spin interaction (of strength A ~ A 2 /<5 a ) associated with 
the Hamiltonian dynamics, whilst mode b mediates the 
collective atomic decay (with rate r b ~ X 2 /n b ). 

The equations of motion for the moments 
{( J x ), (J y ), (J z )}, derived from do not form a 

closed set. However, factorizing the means of operator 
products and taking the limit N — > oo (i.e., neglect- 
ing quantum fluctuations), we obtain a closed set of 
scmiclassical equations, 



X = 2hY -T b ZX, 
Y = -2HX + 2XZX 
Z = -2XXY + T b (X' 2 + Y' Z ), 



T b ZY, 

-2 , v 2\ 



(6a) 
(6b) 
(6c) 



where (X,Y,Z) = ((J x ),(J y ),(J z ))/j with j = N/2, 
and X 2 + Y 2 + Z 2 = 1 (conservation of angular mo- 
mentum). The stable steady-state solutions of (6) ex- 
hibit bifurcations at two critical effective field strengths, 
h c ± = [X± (A 2 - V 2 ) 1,2 ]/2 (we assume A > 0,r 6 ). In 
particular, for h < h°_ and h > h c + the stable steady- 
state solutions are {A ss = Y ss = 0, Z ss = 1}, whereas for 



ht_ < h < h c + one finds 



A ss = ± 



A 2 - 4h 2 
2XA ' 



v - 

1 ss t\ 



Z RR — 
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(7) 
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where A = A + (A 2 -T 2 ) 1 / 2 . Note that at both (supercrit- 
ical pitchfork) bifurcations a detailed stability analysis 
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[T5] shows that a unique steady state becomes unstable 
and two new stable steady states emerge. These semi- 
classical solutions, together with numerical solutions of 
the finite- N master equation ([5]), are plotted in Fig. [2] 
as a function of h/X (note that (J x ) = (J y ) = for the 
finite- TV calculations). The plots indicate both a first- 
and second-order phase transition as a single parameter, 
h, is varied. The first-order (second-order) transition, at 
h = h°_ (h = h c A, involves a discontinuous (continuous) 
bifurcation in X ss and associated behavior in Z ss . Note 
that in the purely Hamiltonian system second-order tran- 
sitions occur at ±/i5_, but the first-order transition has 
no counterpart (for A > 0) and arises here due to a dis- 
sipative instability. The behavior we observe bears some 
relation to critical points found in cooperative resonance 
fluorescence (see, e.g., [IB]). 




FIG. 2: Semiclassical (solid line) and finite- TV steady-state 
inversion and second-order moments for r a /A = 0.01, F b /X = 
0.2, and N = 25 (dotted), 50 (short dash), 100 (long dash). 



To examine this structure and dynamics, we consider 
the transmission of a (weak) probe laser field through the 
medium as a function of the probe frequency, i.e., we ex- 
amine the frequency response of the system. A schematic 
of such a measurement setup is shown in Fig. [TJb). To 
compute the transmission spectrum we retain the two 
cavity modes in our model (i.e., we start from (pj), but 
again perform a linearization (N 3> 1) about the mean- 
field state. We consider the case in which the probe laser 
drives mode 6, and the transmission spectrum T p {y p ) is 
defined as the coherent intensity, at probe frequency v p 
(in the rotating frame), in the output field from mode b. 




FIG. 3: Transmission spectra in the linearised regime, 
for h/X = {-0.6, -0.1, -0.01, ht/X, 0.05, 0.3} (top), and 
{0.5, 0.95, 0.995, h+/X, 1.1, 1.3} (bottom), with microscopic 
parameters K, a /S a = 0.02, X b /X a = 0.32, Kb/5 a = 1, S b — 0, 



giving T a /A = 0.01, Fb/X = 0.05. We set <5 



a,b 



: 0. 



In the large- N limit, quantum fluctuations can be in- 
cluded in the analysis as a first-order correction using a 
large-iV expansion of the Holstein-Primakoff (HP) rep- 
resentation of angular momentum operators [141 . Ap- 
plied in a coordinate system where the mean Bloch vec- 
tor points along the positive z-axis, this takes the form 
J z = N/2 - c+c ~ N/2 and J + = (N - Jcf^c ~ VN c, 
where c (c^) is a bosonic annihilation (creation) operator. 
This linearization about the mean field state leads to a 
master equation of the general form 

p = -i[H m ,,p]+T + D[^]p + T_D[c]p 

+ {T (2cpc-c 2 p- pc 2 ) +H.c.}, (8) 

where H^p is a quadratic in {c, c'~} and the coefficients 
are functions of {h, X, T a , T b } 15J. Eq. Q yields coupled, 
linear equations of motion for (c) and (c^), the eigenval- 
ues of which display a sequence of bifurcations in both 
their real and imaginary parts as h is varied. The phase 
transitions are marked by the real part of one eigenvalue 
going to zero (i.e., critical slowing down) at h = h±. 



In Fig. [3] we plot T p (v p ) (normalized by the maxi- 
mum empty-cavity transmission) for a series of values of 
h around h±. For the chosen parameters, the spectra 
consist of sharp "atomic" resonances superimposed on 
a much broader cavity mode resonance (i.e., K b T b ). 
The locations and widths of the atomic resonances are 
determined by the imaginary and real parts of the above- 
mentioned eigenvalues, respectively. For \h/X\ > 1, the 
main atomic feature is a dip of width 2r& at v ~ 2h, cor- 
responding to a cavity-mediated, collective spontaneous 
emission resonance. For \h/X\ < 1, spin-spin interactions 
play a more significant role and a pair of resonances at 
opposite frequencies feature in the spectrum. As h — * h c + 
both from above and below these two features merge con- 
tinuously into a single peak, centered at v p = 0, which 
ultimately diverges at h = h c + in a pronounced signature 
of the second-order phase transition. The same merg- 
ing and divergence is seen for the first-order transition, 
but only as h — > h c _ from below. For h very small (but 
> lit.), the spectrum consists of two sharp peaks of width 
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~ Tbh/X at v p ~ ±2A. The transition is signaled by a 
discontinuous jump from this two-peaked spectrum to a 
single divergent peak at v p = 0. 

To analyze the entanglement properties of the system, 
we adopt a criterion for bipartite entanglement in col- 
lective spin systems which, for symmetric states, is both 
necessary and sufficient, and reads [TB] 

C v = 1 - (4/A0(Aj2) - (i/N 2 ){J v ) 2 > 0, (9) 

where J v = s'm(ip)J x + cos(tp)J y . Here, we present nu- 
merical results for Cr = max v C v (> 0), which, in fact, 
equals the rescaled concurrence (N — 1)C, where C is the 
two-spin concurrence. In Fig. [4] we plot the steady state 
value of Cr versus h/X, computed from the linearized 
HP model and numerically from ^ for finite N. Both 
transitions are characterized by a sharp peak in the en- 
tanglement at the critical point [5j [6l [7], but they are 
distinguished by a discontinuity in Cr at h = h c _ (for 
N — > oo) as opposed to a discontinuity in dC^/dh at 
h = h c + [T7]. The peaking of Cr at the critical points 
agrees with the conjecture of a general association be- 
tween semiclassical bifurcations and maximal entangle- 
ment in dissipative, non-equilibrium many-body systems 
[18] . In the region where Cr = the state approaches 
a mixture of maximally polarised states possessing large 
fluctuations (see Fig. [2]). Note that, in the adiabatic 
regime considered, the cavity field operator b(t) cx J+{t), 
and so collective spin correlations (and hence Cr) can be 
deduced from moments of the cavity output field, which 
may be measured by broadband homodyne detection. 




- 1 h/X 1 2 

FIG. 4: Maximum entanglement Cr computed from the lin- 
earized HP (N — » oo) model (solid) and from (5) for N = 100 
(dashed), with F a /X = 0.01 and T b /X = 0.2. 

For an experimental realization, we have already men- 
tioned 6 Li in a ring-cavity setup. A suitable system can 
also be designed using the ground states \F = 1, m = ±1) 
of 87 Rb and linearly-polarized cavity modes [TS] . For spe- 
cific parameter values, we consider recent experiments 
with cold atoms inside a high-finesse optical ring cavity 
[19] . i.e., we take gij ~ 2ir ■ 100 kHz and n a ~ 2ir ■ 25 kHz. 
For N ~ 10 6 atoms and a characteristic ratio O^/Aj ~ 
0.0025, we have A Q — In ■ 250 kHz. With a Raman de- 
tuning S a ~ 2ir ■ 2.5 MHz, we then have A ~ 2A„/<5 a ~ 
2tt-25 kHz and T a ~ K a (X a /S a ) 2 ~ 2tt-0.25 kHz. Ground 
state magnetic level shifts of tens of MHz would suffice 
to ensure distinct Raman channels. Mode b may be more 
strongly damped (i.e., the two cavity polarizations have 
different finesses), e.g., Kb ~ 2ir ■ 250 kHz, and, with 



Xb — 2ir ■ 25 kHz and 8b — 0, we would then have Tb — 
X 2 / Kb — 2ir ■ 2.5 kHz 3> r a . Finally, the rate for single- 
atom spontaneous emission (neglected in our model) is 
estimated by r at nf 3 ./(4A?) < 2ir ■ 0.01 kHz < A,r fc for 
an atomic exited state decay rate T at = 2ir ■ 6 MHz. 

To conclude, we have proposed a feasible cavity QED 
system that is described by a dissipative LMG model 
and exhibits both first- and second-order non-equilibrium 
quantum phase transitions as a function of a single effec- 
tive field parameter. Measurements on the cavity out- 
put light fields provide quantitative probes of the critical 
behavior. The system also offers opportunities for in- 
vestigating phase transitions in response to variation of 
the strength of dissipation (i.e., Tb), for studying time- 
dependent behavior, such as entanglement dynamics, and 
for preparing very highly entangled states, which typi- 
cally occur for short interaction times |15j and may in 
principle be "frozen" by switching off all optical fields. 
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